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We study non-adiabatic charge pumping through single-level quantum dots taking into account 
Coulomb interactions. We show how a truncated set of equations of motion can be propagated 
in time by means of an auxiliary-mode expansion. This formalism is capable of treating the time- 
dependent electronic transport for arbitrary driving parameters. We verify that the proposed method 
describes very precisely the well-known limit of adiabatic pumping through quantum dots without 
Coulomb interactions. As an example we discuss pumping driven by short voltage pulses for various 
interaction strengths. Such finite pulses are particular suited to investigate transient non-adiabatic 
effects, which may be also important for periodic drivings, where they are much more difficult to 
reveal. 
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I. INTRODUCTION 

In 1983 Thouless [I] proposed a simple pumping mech- 
anism to produce, even in the absence of an external bias, 
a quantized electron current through a quantum conduc- 
tor by an appropriate time-dependent variation of the 
system parameters. Experimental realizations of quan- 
tum pumps using quantum dots (QDs) were already re- 
ported in the early 90's [HE]. More recently, due to the 
technological advances in nano-lithography and control, 
such experiments have risen to a much higher sophistica- 
tion level, making it possible to pump electron [3H5] and 
spin [7] currents through open nanoscale conductors, as 
well as through single and double QDs [811 . 

Early theoretical investigations where devoted to the 
adiabatic pumping regime within the single-particle ap- 
proximation jl2H14] . This is well justified for experi- 
ments with open QDs, where interaction effects are be- 
lieved to be weak 115] and the typical pumping param- 
eters are slow with respect the characteristic transport 
time-scales, such as the electron dwell time t^- This time- 
scale separation enormously simplifies the analysis of the 
two-time evolution of the system. Within the adiabatic 
regime, inelastic and dissipation |16| effects of currents 
generated by quantum pumps were analyzed. Further- 
more, issues like counting statistics [T7], memory effects 
[18], and generalizations of charge pumping to adiabatic 
quantum spin pumps were also proposed and studied [151 - 
[II]. 
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Non-adiabatic pumping has been theoretically inves- 
tigated within the single-particle picture, either by us- 
ing Keldysh non-equilibrium Green's functions (NEGF) 
with an optimal parametrization of the carrier operators 
inspired by bosonization studies 22 , or by a Flouquet 
analysis of the S'-matrix obtained from the scattering ap- 
proach [23i . While the first approach renders complicated 
integro-differential equations for the Green's functions as- 
sociated to the transport, the second one gives a set of 
coupled equations for the Flouquet operator. It is worth 
to stress that, in both cases the single-particle picture is 
crucial to make the solution possible and it is well estab- 
lished that both methods are equivalent [Ml US] . 

Several works have provided a quite satisfactory de- 
scription of quantum pumping for weakly interacting sys- 
tems. In contrast, the picture is not so clear for situa- 
tions where interaction effects are important. Different 
approximation schemes have been proposed to deal with 
pumping in the presence of interactions and to address 
charging effects, which are not accounted for in a mean- 
field approximation. Typically, two limiting regimes have 
been studied, namely, the one of small pumping frequen- 
cies ri, such that r^Td ^ 1 (adiabatic limit) [2SH5^ and 
the one of very high frequencies, fird ^ 1 (sudden or 
diabatic limit) [55HS5] . Nonadiabatic pumping is mainly 
studied as a side effect of photon-assisted tunneling [36l - 
155] . where Q.Ti » 1. 

Unfortunately, it is quite cumbersome to calculate cor- 
rections to these limit cases. For instance, the analysis 
of higher-order corrections to the adiabatic approxima- 
tion for the current gives neither simple nor insightful 
expressions [32]. In addition to the theoretical interest, a 
comprehensive approach bridging the limits of fird ^ 1 
and f2Td ^ 1 has also a strong experimental motiva- 
tion: Most current experimental realizations of quantum 
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pumping deal with QDs in the Coulomb blockade regime 
and r^Td ~ 1. This regime was recently approached (from 
below) by means of a diagrammatic real-time transport 
theory with a summation to all orders in fl '39]. However, 
the derivation implied the weak tunnel coupling limit, 
whereas experiments [40ti44] typically rely on tunnel cou- 
pling variations which include both weak and strong cou- 
pling. 

To address the above mentioned issues and to ac- 
count for the different time scales involved it is natural 
to use a propagation method in the time domain [^51 - 
I49j . In this work we express the current operator in 
terms of density matrices in the Heisenberg representa- 
tion. We obtain the pumped current by truncating the 
resulting equations-of-motion for the many-body prob- 
lem. The time-dependence is treated exactly by means 
of an auxiliary-mode expansion [37| \SU\ . This approach 
provides a quite amenable path to circumvent the usual 
difRculties of dealing with two-time Green's functions 
[50] . Moreover, it has been successfully applied to sys- 
tems coupled to bosonic reservoirs [51] and to the descrip- 
tion of time-dependent electron-transport using general- 
ized quantum master equations for the reduced density 
matrix |47 1 I52 [ [53]. Since the auxiliary-mode expansion is 
well controlled [SI] , the accuracy of our method is deter- 
mined solely by the level of approximation used to treat 
the many-body problem. 

The formalism we put forward is illustrated by the 
study of the charge pumped through a QD in the 
Coulomb-blockade regime by varying its resonance en- 
ergy and couplings to the leads. The external drive is 
parametrized by a single pulse, whose duration and am- 
plitude can be arbitrarily varied. By doing so, the for- 
malism is capable to reproduce all known results of the 
adiabatic limit and to explore transient effects beyond 
this simple limit. 

The paper is organized as follows. In Sec.[ll]we present 
the resonant-level model, as well the theoretical frame- 
work employed in our analysis. In Sec. [Illj we introduce 
the general propagation scheme, suitable to calculate the 
pumping current at the adiabatic regime and beyond 
it. Next, in Sec. |IV[ we discuss few applications of the 
method. Finally, in Sec. [V] we present our conclusions. 



II. TIME-DEPENDENT INTERACTING 
RESONANT-LEVEL MODEL 



resonance model. In addition, in the Coulomb blockade 
regime F is much smaller than the resonance charging 
energy U. 

A. Setup 

The total Hamiltonian is given by the usual threefold 
decomposition into a quantum dot Hamiltonian i?dot, a 
Hamiltonian iJicads representing the leads, and a coupling 
term i/coup, namely 

H ~ Hdot + Hicads + -ffcoup ■ (la) 

The QD is modeled by a single level of energy ed{t), 
which can be occupied by spin-up and spin-down elec- 
trons, which interact through a contact interaction of 
strength U. The QD Hamiltonian reads 

-ffdot = ^ £d{t)ns + Ufi^fii , (lb) 

where Us = clcg, and Cg stand for electron number, 
creation and annihilation operator for the respective spin 
state s =t, I in the dot. 

The two reservoirs, labeled as L (left) and R (right), 
are populated by non-interacting electrons, whose Hamil- 
tonian reads 

i^lcads = Y^c,k{t)Vaks^aks , (Ic) 

Q(iL,R ks 

where {b^^f^g\ and {baks\ stand for the electron creation 
and annihilation operators for the a-reservoir state A:s, 
respectively. The reservoir single-particle energies have 
the general form e^k {t) — £qA; + [t) with the Aeq, ac- 
counting for a time-dependent bias. The stationary cur- 
rent due to a time-dependent bias was already addressed 
several years ago [55]. For pumping, we take l^ea{t) = 0, 
as usual. 

Finally, the coupling Hamiltonian is given by 

^coup = E E ^^l^A + H.C. , (Id) 

ak s 

with {TqA;} denoting the coupling matrix element be- 
tween the QD and the reservoir a. 



The standard model to address electron trans- 
port through QDs is the Anderson interacting single- 
resonance model coupled to two reservoirs, one acting 
as a source and the other as a drain. Despite its simplic- 
ity, the model provides a good description for Coulomb- 
blockade QDs and for QDs at the Kondo regime, where 
the electrons are strongly correlated. In this paper we 
address the Coulomb-blockade regime, for QDs whose 
typical line width F is much smaller than the QD mean 
level spacing 5, justifying the use of the Anderson single- 



B. Equation of motion approach 

We are interested in the electronic current from reser- 
voir a to the QD state s, which can be obtained from 
current operator 



Ls{t) = iY,[T^k{t)bikAi)^s{t) 



(2) 
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Here and in the following we use units where the ele- 
mentary charge e = 1 and the reduced Planck constant 
h = 1, unless otherwise indicated. To calculate Jas{t) 
we use the following equations of motion, which are ob- 
tained from the Hamiltonian [Eqs. ([I])] by means of the 
Heisenberg equation, 



idtcsit) = es{t)cs{t) -f Ucs{t)hs{t) 



ak 



ak 



+nk{t)c\{t)hc.ks{t) 



(3a) 
(3b) 

(3c) 



Analogous equations hold for V^f.^ and c\. 

In the spirit of the scheme introduced by Caroli and co- 
workers [51) , we assume an initially uncorrelated density 
operator of the combined system, i. e., we set Takito) 
for to —oo. Further, we apply the so-called wide-band 
limit [57 , where the square of the tunneling element Tak 
is inversely proportional to the density of states pa at 
energy Sak ■ By means of the lead Green function [5 7) 



9ak{t,t') = exp 



dt"e„k{t") 



we can define the decay rate 



(4) 



(5a) 



which becomes local in time in the wide-band limit, 
namely 

ra(i,t')= j dekPiek)T^k{t)9c.k{t,t')T*k{t') 

^T„{t)S{t-t') . (5b) 



In the following we replace the sum in Eq. (5a I by the 
expression involving the (5-function in Eq. (5b I. 



The equation of motion for the reservoir operators 



baks, Eq. (3a), is now readily integrated, yielding 

t 

kks{t)^B^ksit)-i J dt' g^k{t.t')T^k{t')cs{t') , (6a) 

to 

where we have used the lead Green functions, Eq. ([4]), 
and introduced 



Baks{t) = gak{t,to)baks{to) ■ 

Equations ([6| are used to rewrite Eq. ( 3b ) as 



(6b) 



idtCsit) 



-J2nk{t)Baks{t). 
ak 



Csit) 



(7) 



Here the wide-band limit, Eq. ([5]), is employed to ob- 
tain the decay term, proportional to T{t) = '^a^ait). 
Similarly, we can rewrite Eq. (|3c| as 



(8) 



idtn-s(t) = ^ [-T^k{t)Bl,,^{t)c-s(t) 

ak 



+T:^{t)cUt)B^ks{t)\ -ir(i)n,(i) . 

Here again the time integral of baks{t) is reduced to a 
decay width due to the wide-band limit 

C. Expectation values and truncation schemes 

The expression for the time-dependent current is given 
by the expectation value of the current operator J^g de- 
fined in Eq. ([2| . As will become clear later on, it is useful 
to write this expectation value as 

•^"W = i E (-^-w) = y En-W (9) 

s s 

with the current matrices of the first order 

n^S)^iY.^ck{t)(bl,^{t)cs{t)) . (10) 

k 

These current matrices are an essential ingredient of our 
propagation scheme, which is based on finding equations 
of motion for Has- Such equations have been derived 
starting from a NEGF formalism for non-interacting elec- 
trons [50] . 

Exactly as for the operator equations above we can use 
baks from Eq. ^ and employ the wide-band limit ([5| for 
the current matrices defined in Eq. (10 1. This leads to 
the following decomposition 

Uasit) =Ksit) + J2^o.k{t)Kksit) , (lla) 

k 

Ksit) = (nsit)) , (lib) 

Kks{t)=UBlkJt)cs{t)) . (11c) 



Having derived all relevant equations of motion for the 
operators we can specify the respective equations for the 
two contributions H' and H". The term H' is the simplest 
and is basically given by the equation of motion for hs , cf. 
Eq. ([8|. The corresponding equation for the occupation 
ns{t) = {hs{t)) reads 

dtusit) = 2Re Y,Tc.k{t)Kksit) - r(i)n,(t). (12) 



The above relation can be viewed as the charge conserva- 
tion equation for the QD. The rate by which the charge 
in the QD changes is equal to the total electronic cur- 
rents. The first term at the r.h.s. of the equation can be 
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interpreted as the current flowing into the QD, whereas 
the second term gives the current flowing out. 

Since we do not consider a spin-dependent driving or 
spin-polarized initial states it is ns{t) — ns{t). This re- 
lation is not explicitly used in the derivation, but is em- 
ployed as a consistency check throughout the analysis. 

The evaluation of H" requires the solutions for both, 
the lead operator baks and the dot operator Cg. Using 
those, we write 



Here we have introduced the abbreviation 



Aaks{t) = e°fe 



(13) 



(14) 



and used that 



fa{£k) = fak 



(15) 



with fa (e) the Fermi function describing the equilibrium 
occupation of lead a. The last term in Eq. (13) uses the 



auxiliary current matrices of the second order 
^aks^i(Bl^At)csit)n-sit)] 



(16) 



which will be subject to further approximations in the 
following. 

Before we turn to the approximations, we would like 
to briefly discuss the physical meaning of ^aks- The 
equation of motion for the two-electron density matrix 
{ns{t)hs{t)) reads 

dt {h,it)hs{t)) = - 2T{t) {hsit)hs{t)) 

+ 2Re Y,Tc.ks{t)^cks{t) , (17) 

aks 

which follows from Eq. ([s]). The two-electron density 
matrix may be interpreted as the occupation of one 
quantum-dot level under the condition that the other one 
is occupied. The rate of change of this conditional occu- 
pation is consequently given by tunneling into and out of 



the respective dot state under the same condition. The 
latter process is described by the first term on the r.h.s. 
of Eq. (17). The former process is governed by the aux- 
iliary current matrices ^aks, which can be rewritten in 
the suggestive form 



2Re ^T„fc,(t)$„fe,(t) = /.Ls(.t)hrs{t) 



(18) 



Consequently, the current matrices <&afcs describe the 
conditional current from reservoir a into the quantum- 
dot level with spin s. 

1. Hartree-Fock approximation 

The simplest approximation to ^aks consists in using 
the following factorization 



^^lit)^i{Bl,^{t)cAt)){n,{t)) 
= n-s{t) Kksit) 



(19) 



Inserting this expression into Eq. (13), results in the fol- 
lowing equation of motion 

dtK^M = i [Aofc.(t) - un-4t)] n::,,(t) 

+ T*kit)fc.k. (20) 

This result is equivalent to the Hartree-Fock approxi- 
mation applied to the Anderson model standard two- 
electron Green function [57| . As any mean field approach, 
it does not lead to a double resonance Green function, 
which is required to properly account for charging ef- 
fects. Hence, as it is well known, a good description of 
the Coulomb-blockade regime requires going beyond this 
level of truncation in the equations of motion. 



2. Hubbard I approximation 

Instead of factorizing ^aks directly, we proceed by de- 
riving its equation of motion. By means of Eqs. ^ we 
get 



dt^ckAt) = i [Aafcs ~U- iVit)] $„fe,(t) + J2 T^'k'it) {Blksit)Ba'k'sit)nrsit) 

a' .k' 

+ J2 [^"''^■'(*) {BLsit)^sit)Bl,,,m{t)) - r:,,,(i) (^Bi,At)cs{t)cUt)Ba,k',{t) 



(21) 



a'k' 



Note that the term proportional to U has only four op- The approximation consists in neglecting matrix ele- 

erators in the expectation values because of figfig — fig. ments involving opposite spins, which renders the follow- 
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ing factorizations 



baksi'^o)ba' k' sito)ns{t)) ~ fakSaa'Skk'Tlsit) , 



blksiio)cs{t)bl^,,{h)c-s{t) 
blksito%it)cUt)b^>k's{k) 



(22) 



This approximation for the density matrices is equiva- 
lent to the truncation scheme employed in the NEGF 
approach used for the study of Coulomb blockade regime 
(high-temperature limit of the Anderson model) [SJ • 

As a result of the factorization, we obtain the follow- 
ing compact equation of motion for the approximated 
second-order current matrices 

dt^aksit) = i [A^ks ~U~ iT{t)] l>„fc,(0 

+ T*,it)Uknsit) . (23) 

The equations of motion for ns{t) [Eq. ([T2|], Il'^ks(t) [Eq- 
( fl3t with ^aks replaced by <^aks] and for ^aks{t) [Eq. 



(23l] form a closed set of equations, which can be solved 



by means of an auxiliary-mode expansion discussed be- 
low. 



III. AUXILIARY MODE PROPAGATION 
SCHEME 

The general idea of the auxiliary-mode expansion con- 
sists in making use of a contour integration and the 
residue theorem to perform the energy integration, for 
instance, in Eq. (11). To this end the Fermi function is 



expanded in a sum over simple poles (or auxiliary modes) 
and the respective integrals are given as finite sums, cf. 
Appendix |A] 

The transition to auxiliary modes (denoted by the in- 
dex p) is facilitated by the following set of rules 



E 



£ak{t) 
'^akit)fak 

T^k{t)Kksit) 



Xapi^) 1 

1 



(24a) 
(24b) 



r„(i)+T„(i) 



Vn 

/ J asp 



(24c) 

which are derived in Appendix]^ The first rule replaces 
the reservoir energy Sg k by the (complex) pole xtp of 
the expansion, cf. Eq. (A5). The second rule replaces 



the Fermi function by the respective weight, which is 
the same for all auxiliary modes. Finally, the third rule 
provides the actual expansion for the current matrices. 
Applying these rules, the current matrices become 



Ksit) 
KM 



Tait) 



2 

r„(t) 



-nsit), 
Tgit)J2Kspit)- 



(25a) 
(25b) 

(25c) 



The equation of motion for the auxiliary matrix n"^,^ is 
obtained from Eq. (13). One arrives at 

idtKspit) = [esit)-im/2-xlpit)]Kspit) 

+ ^Tgit) + U ^a^spit) . (26) 

The equations of motion for the auxiliary matrices ^asp 
are quite similar to those of Eq. 



(25c 



idt<^>aspit) 



sit) + U-xlJt)- 



namely, 

3r{t) 



To.it) ns{t) 



^aspit) 

(27) 



The solution of the above equations still requires a com- 
plete description of the population dynamics given by 
ns{t). The latter can be directly obtained from Eq. (12) 
in terms of the current matrices 



dtus it) = -r(t) Us {t) + 2Re Ks it) 



(28) 



This concludes the derivation of the auxiliary mode prop- 
agation scheme. The set of equations ( p6| to (28 1, 
with initial conditions ns(io) = 0, n^j,p(tj^= 0, and 

^aspito) — 0, can be solved numerically using standard 
algorithms. Before the desired time dependence of the 
parameters Ssit) and V{t) sets in, the system has to 
be propagated until a steady state is reached. In this 
way, transient effects arising from the choice of the initial 
state are avoided. For convenience we derive in Appendix 
|B]the expressions for the stationary occupations, which 
may also be used as initial values for n^. 



IV. NON-ADIABATIC PUMPING 

In this section we present two applications of the for- 
malism developed above. As shown below, one of the 
interesting features of non-adiabatic pumping is an in- 
creasing delay in the current response to the external 
drive with growing driving speed. Hence, in distinction 
to the adiabatic limit, the current caused by a train of 
pulses can show interesting transient effects, whenever 
the pulse period is shorter than the system response time. 
To better understand non-adiabatic driving effects, we fo- 
cus our analysis on single pulses and vary the speed by 
which their shape is changed. 

It is worth stressing that our propagation method does 
not possess restrictions on the time dependence of the 
system driving parameters. In other words, the external 
time-dependent drive can be just a single pulse or a train 
of pulses, it can also be either fast or slow as compared 
with the system internal time scales. 



A. Symmetric monoparametric pumping 

Let us begin by discussing the current generated by 
a single Gaussian voltage pulse changing the resonance 
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energy as 

edit) = So + £iexp [-(t/tp)^] 



(29) 



p sets the pumping time-scale. We take Tait) to 



Here . . _ 

be time-independent and equal for both leads, Ti^{t) — 
^nit) = To/S. Since thereby Jl = Jr, we will consider 
only Jl in the following. 

Figure [T^ shows the time dependence of the resonance 
energy according to Eq. (29). The two bottom panels 
show the instantaneous current Jl as a function of time 
for both the non-interacting (U = 0) and the interacting 
{U ^ 0) case. In the hmit of large tp, we use as a check 
for our results an analytical expression for the pumped 
current Jl, obtained for U = within the adiabatic ap- 
proximation [T^ I28| . 

Here, due to the L/R symmetry, there is no net charge 
flowing through the QD. At any given time both leads 
pump the same amount of charge in or out. In the driv- 
ing scheme defined by Eq. ( 29 ) , the QD is initially nearly 



empty. At t — 0, the resonance energy favors an almost 
full occupation. For very slow pumping, large tpFo, the 
current Jl depends only on the resonance energy ed{t): 
As the resonance dives into the Fermi sea, the QD is 
loaded with charge and the process is reversed as ed{t) 
starts increasing. This is no longer true when the drive 
is faster and tpTo decreases: Now one observes a retar- 
dation effect, namely, the Jl depends not only on the 
resonance position, but also on driving speed. For fast 
driving one needs to integrate Jl over times much longer 
than tp to observe a vanishing net charge per pulse. 



B. Constrained two-parameter pumping 

The pumped currents JL.R(t) characterize the time- 
dependent electron response to the external drive. How- 
ever, in most applications one is only interested in the 
charge pumped per cycle Qc or per pulse Qp. In the 
latter case, Qp is given as time integral over the current 
which we write in a symmetric way 



+00 



Mt)] 



(30) 



One of the beautiful lessons learned form the investi- 
gation of adiabatic pumping, establishes a proportion- 
ality relation between Qp and the area swapped by the 
time-dependent driving forces in parameter space [12j . In 
other words, the total charge flowing through a QD per 
cycle (or per pulse) in a single-parameter adiabatic pump 
vanishes. Due to the constraints of single-parameter 
pumps, in most applications at least two parameters are 
used [2HS1 [5HTH HOPIj. On the other hand, by using 
a single-gate modulation one can realize a constrained 
two-parameter pump |43j . which implies that the time- 
dependence of the parameters is ultimately coupled due 
to the modulation of only a single gate voltage. In the 
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FIG. 1: a) Time dependence of the resonance energy accord- 
ing to Eq. ( 29 1 ; b) and c) Time-resolved current Jl for dif- 



ferent pulse lengths tp and U = and U/To ~ 10. Param- 
eters used: eo = To, ei = — 2ro, fia ~ 0, ksT = O.lFo and 
A^F = 160 (number of auxiliary modes). Dots denote the 
adiabatic limit 28 . 



following we will investigate the implications of this sce- 
nario for non-adiabatic pumping. 



1. Pulse Scheme 
Specifically, let us consider voltage pulses of the form 
Sit,5) = l-2exp\-{2Vh^t/tp-S)'^] . (31) 
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Here tp measures the characteristic pulse time, whereas 5 
governs the time the pulse sets in. The numerical factor 
2y/\n 2 ensures that tp is the full width at half maximum 
of the pulse, which simplifies the following discussion. By 
tuning the delay one can conveniently switch between a 
single parameter [S = 0) and a two parameter setup {6 ^ 
0). Further, the time-dependence of the resonance energy 
and the coupling strengths (decay widths) are chosen as 



Edit) = eo + eiS{t,0) , 

ro 

2 
2 



rR(t) 



l-5(t,0)] , 



(32a) 
(32b) 

(32c) 



This choice takes into account that the coupling strengths 
depend exponentially on the gate voltage [321 ■ The con- 
straint is imposed by setting (5r = and the specific value 
of For this driving parameterization, the resonance 
and the decay widths are Eq and Fq, respectively, for both 
asymptotic limits of |t| ^ tp. In the following, the pa- 
rameters are taken as Eq — 0, Ei/Tq — 4, kBT/To = 1/10, 
and interaction energy either U — or U/Tq = 10. In 
Fig. [21 the time dependence of Ed and Fl/j^ is illustrated 
for tnree cases (5l = and — ±1. As mentioned 
above, in each case the coupling to the right reservoir 
Fr follows the time dependence of the resonance energy. 
When the latter attains its minimal value at t ^ 0, which 
brings the energy well below the chemical potential of the 
reservoirs, the coupling to the right reservoir is minimal. 
On the other hand, the behavior of the coupling to the 
left reservoir can be influenced by the value of (Sl- For 
Si^ — —1 the maximum of Fl comes before i = 0, while 
for ~ +1 it is attained after t = 0. In the case (5l = 
the coupling to the left reservoir is maximal simultane- 
ously with Fr being minimal a,t t — 0. In the following 
the response to these drivings will be investigated. 
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FIG. 2: Time-dependence of the resonance energy Ed (upper 
row) and the decay widths Fl (lower row, blue/full line) and 
Fr (lower row, red/broken line) for three different cases: a) 
(5l = — 1, b) 5l = 0, and c) 5l = 1. The dotted lines indi- 
cate the chemical potential in the reservoirs and the shaded 
area shows the times when the resonance energy is below the 
chemical potential. 



2. Adiahatic Pumping 

Knowing the time dependence shown in Fig. [2] one can 
readily predict the behavior of Qp in the adiabatic limit. 
In this case, electron flow occurs when the resonance en- 
ergy matches the chemical potential of the reservoirs. In 
our pulse scheme, ed(^) equals the chemical potential at 
t = —tp/2 and t = +tp/2 corresponding to the onset of 
charging and de-charging of the QD, respectively. Fur- 
ther, the direction of the net current is determined by the 
difference of the couplings to the reservoirs at these very 
times. For example, for (5l = — 1 one finds Fl > Fr while 
charging and Fl < Fr while dc-charging. Consequently, 
the net current is directed from left to right and Qp is 
expected to be positive. For (5l = +1 the situation is 
opposite and Qp should be negative. Finally, for ^l = 
the couplings are equal at both instants of time and the 
net current is vanishing. These expectations are con- 
firmed by our results for the adiabatic regime, tpFp S> 1, 
and different values of CZ/Fq, which are shown in Fig. [3^. 
As already mentioned, one observes Qp = for (5l = 
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(monoparametric pumping). As begins to increase, 
IQpl increases as well. In this scenario, when the reso- 
nance energy matches the chemical potential, electrons 
load the dot from the left (or right) and later they are 
unloaded to the right (or left, depending on the sign on 
S). For larger values of |(5l|, the left reservoir participates 
less in the loading or unloading of the QD and the charge 
per pulse vanishes accordingly. 

For interaction strengths C/ ^ ei the double occupa- 
tion of the QD is suppressed and consequently, in the 
adiabatic regime, Qp is half the value of Qp for the non- 
interacting case. The numerical results indicate that 
within the Hubbard 1 approximation, [/ ^ does not 
introduce new time scales to the problem for tpFo ^ 1, 
and its major effect is to correct the spin degeneracy fac- 
tor in the equations for the C/ = case. 



3. Non-adiabatic Pumping 

None of the aforementioned features are observed in 
the non-adiabatic pumping regime. Figure [3]d shows, for 
example, that, for short pulses there is no simple relation 
between Qp for U = and for U 0. Moreover, com- 
pared to the adiabatic regime the charge per pulse can 
be substantially larger in this regime. Unfortunately, the 
behavior of Qp in this regime is not as easily predicted 
in general, since the evolution of the parameters {ed, Fl, 
Fr} after the onset of loading and unloading has to be 
taken into account. This is because in the non-adiabatic 
regime the QD charging and de-charging is delayed with 
respect to the external system changes, as it was shown 
in Sec. |1V A[ Taking, for example, the case (5l = one 
finds from Fig. [2]d, that Fl > Fr while the resonance 
energy is below the chemical potential and charging oc- 
curs. During the de-charging, when Ed > /^l,r, one finds 
Fr > Fl. Consequently, the current is expected to flow 
mainly from left to right, which leads to a positive charge 
per cycle. This is confirmed by the results shown in Fig. 



|3]d. The quantitative behavior depends on the precise 
magnitude of the delay, which is determined by the pulse 
length. However, from the analysis presented above and 
for sufficiently short pulses one concludes that Qp has to 
be positive independent of ^l- The interesting implica- 
tions of this result will be discussed at the end of this 
section. 

Finally, in Fig. |4] we summarize and corroborate the 
discussion of the non-adiabatic pumping. It shows the 
charge pumped due the pulse as a function of pulse length 
tpTo in the non-interacting {U = 0) and the Coulomb 
blockade regime {U = IOFq). In the latter case Qp < 1 
for all pulse lengths. As discussed above, the amount of 
pumped charge Qp depends very strongly on the value of 
ipFo. In the limit of large pulse lengths, Qp approaches 
the respective adiabatic value, while for tpFo — >■ the 
charge per pulse vanishes. Moreover, one finds that Qp 
is indeed positive for small pulse lengths. This has the 
intriguing consequence that the charge per pulse can 
change its sign sweeping from short to long pulses. This 
is shown in Fig. |4j3 for i5l = 1, where Qp is negative in 
the adiabatic regime. A more general and quantitative 
analysis of this effect is certainly desirable, but beyond 
the scope of this article. It may lead, however, to in- 
teresting new applications. It is also worth to mention, 
that by changing the pumping parameters it is possible 
to optimize the charge pumped per pulse and in partic- 
ular to find situations where Qp = 1, which may be very 
interesting for metrology purposes |58| . 



V. CONCLUSIONS 

We presented a new method to analyze non-adiabatic 
charge pumping through single-level quantum dots that 
takes into account Coulomb interactions. The method is 
based on calculating the time evolution for single-electron 
density matrices. The many-body aspects of the problem 
are approximated by truncating the equations of motion 
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one order beyond mean field. The novelty is the way the 
time evolution is treated: By means of an auxiliary-mode 
expansion, we obtain a propagation scheme that allows 
for dealing with arbitrary driving parameters, fast and 
slow. The method presented in this paper can be applied 
to a wide range of coupling parameters , provided one 
avoids the Kondo regime. Hence, we are not restricted 
to the weak coupling limit where Qp, the charge pumped 
per pulse, is rather small. 

The presented results for single-pulses are also valid 
for pulse trains, provided the time between the pulses 
is sufficiently long. One can expect to find qualitatively 
new and interesting effects by decreasing the time lag. 
The propagation scheme allows, in principle, for study- 
ing transient effects. In addition, by propagating over a 
periodic sequence of pulses it constitutes a complemen- 
tary approach to the more familiar periodic driving. In 
this regard, our propagation scheme has the potential to 
be a valuable tool and provide deeper insights into non- 
adiabatic quantum pumps. 
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Appendix A: Auxiliary- mode expansion 



Here we motivate the rules given in Sect. |ni| To begin 
with we introduce correlation functions, which can be 
approximated by finite sums. Then we write the current 
matrices in terms of these finite sums. 



1. Correlation functions and mode expansion 

As we will show later, in the present case we have to 
consider the following reservoir correlation function 

C^{t\t) = ^T„fc(t)g„fc(t',t)T^,(i')/a(£fe) 



de 
2^ 



expi i J dt"[e + Aea{t")] 



(Al) 



where the line-width function Ta is defined as usual [57j 
r„(£, t',t) = 2nY, Tak{t)T*f,{t')6{e - Sak) 

k 

= 2nT^it)T*{t')p^. (A2) 
In the second line we have used the wide-band limit. 



In order to perform the energy integration in Eq. ( Al ) 



we expand the Fermi function /(e) as a finite sum over 
simple poles 

(A3) 



2 P±-?V£-Xp £-Xp 



with = fJ-ixpf/S and ImXp > 0. Instead of using the 
Matsubara expansion [59], with poles Xp = i7r(2p— 1), we 
use a partial fraction decomposition of the Fermi func- 
tion [54 , which converges much faster than the standard 
Matsubara expansion. In this case the poles Xp = -iilJl^ 
are given by the eigenvalues Zp of a A^f x matrix [54] . 
The poles are arranged such that all poles Xp (Xp ) are in 
the upper (lower) complex plane. As in the Matsubara 
expansion all poles have the same weight. 



Employing the expansion given by Eq. ( A3 1 , one can 



evaluate the energy integrals by contour integration in 
the upper or lower complex plane depending on the sign 
of t — t' . Thereby, the integral in Eq. (All becomes a 



(finite) sum of the residues. For t>t' one gets 

c^{t',t) = h:^{t)6{t-t') 

+ T^{t)Y,Cap{t,h) , (A4a) 

(A4b) 

with the auxiliary modes for reservoir a given by 



xipW = K + A£„(i)]+ 



(A5) 



Here, /Xq, is the chemical potential and Aea{t) is due to 
the time-dependent single-particle energies Sak (t) of the 
reservoir Hamiltonian [Eq. (Ic)]. 



Current matrices 



The set of equations (13) and (23) can be formally 



solved. In order to write down these solutions we define 
the following functions 



(A6a) 
(A6b) 



With these definitions the formal solution of Eq. ( 13 ) 
reads 



Jto 

X T*,{t')Uek)^iU^o.ks{t')] , (A7) 

where we have assumed ^^ksi^o) = Oj corresponding to 
our choic e of an initially uncorrelated density matrix (see 
Sec. II A). An analogous equation holds for ^afcs(i), again 
with $afcs(to) = 0. We can combine these two expres- 
sions to get for the second part of the current matrix 
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z z z 

KAt) ^T.'^o.k{t)KkAt) = J dt'C^{t',t)gsit,t') -iuj dt'gs{t,t') J dt"C^{t",t)9]!{t'X)n,it") , (A8) 



wher e we have used the definition of the correlation function Cq given by Eq. (All. Finally, by means of the expansion 

(A9a) 



( A4 ) of the correlation functions we obtain an expansion of the current matrices 



t t t' 

Kspit) = J dt'Co.p{t',t)g,{t,t')^iU J dt'gs{t,t') J dt"Cc.p{t" ,t)g^ {f ,t")ns{t") , (A9b) 

to to to 



which resembles the last rule of Eqs. (24 1. Using the ex- Equation (B3) is a non- linear equation for rig and has to 



plicit expression for II^^p and taking the time derivative 
one can easily verify the first two rules given by Eqs. ( 24 1 . 



Similarly, one also obtains an expression for ^asp, which 
reads 



EE / dt"Capit",t)g^it,t")n,{t") . (AlO) 



The time derivative of this expression is given by Eq. 
(127). 



Appendix B: Stationary occupations 

If neither the couplings Tak (and thus F) nor the levels 
£s or Eak depend on time the level occupations Us and 
the currents Jq, converge to stationary values. These val- 
ues can be obtained by setting all time derivatives in the 
respective equations of motion to zero. In order to sim- 
plify the notation we characterize the stationary values 
by omitting the time argument. 

We will specify below the level occupations for the two 



approximations discussed in Sec. II C Therefore we use 
= ^2Re^T„,n::,,, (Bl) 



ak 



which follows directly from Eq. ([TT[3). 



Hartree-Fock 



Within the Hartree-Fock approximation [Sec. II C 1 
get from Eq. ( 20 ) 



n" 



A 



aks 



(B2) 



Plugging this into Eq. (Bl), changing the k summation 



into an integral over e and using the definition (14), we 
get for the wide-band limit [Eq. ^] 

de fa{e) 



2n 



(e-es-Unrsf + {^y 



(B3) 



be solved numerically. 



2. Hubbard I 

We obtain the stationary conditional current ^aks for 
the Hubbard I approximation [Sec. II C 2 from Eq. ( 23 ) 

as 



T^kfakng 



;7 + iF 



(B4) 



This expression can be used for the stationary H" in Eq. 
p) 



tt// _ • ^Oik J Q^' 
aks A 

^aks 



T*k fak Ung 



^aks [Aafcs — U + IF] 



(B5) 



We use Eq. ( |B1[ ) and the definition ( |14[ ) and finally get 

iwing integral 

Us) [A'{e)+n,A"{e)] 



for the occupation the following integral 
de 



27r' 



A'{s) 

A"{e) ^ A'{e) 



1 



U [4(e-£,) - U] 



(B6a) 
(B6b) 

(B6c) 



This time the equation is linear in rig and can be solved 
explicitly. In the limits [/ — and U — oo it is 
A"{e) = and A"{e) = ~A'{e), respectively. The for- 
mer limit corresponds to non-interacting electrons and 
Eq. ( |B6a ) gives the correct expression for the occupation 
|57) . The latter case describes the situation with very 
strong interactions. 
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